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ABSTRACT 

We present a systematic study of line widths in the [0 III] A5007 and Ha lines for a sample of 
86 planetary nebulae in the Milky Way bulge based upon spectroscopy obtained at the Observa- 
torio Astronomico Nacional in the Sierra San Pedro Martir (OAN-SPM) using the Manchester 
Echelle Spectrograph. The planetary nebulae were selected with the intention of simulating 
samples of bright extragalactic planetary nebulae. We separate the planetary nebulae into two 
samples containing cooler and hotter central stars, defined by the absence or presence, respec- 
tively, of the He II A6560 line in the Ha spectra. This division separates samples of younger and 
more evolved planetary nebulae. The sample of planetary nebulae with hotter central stars has 
systematically larger line widths, larger radii, lower electron densities, and lower H/3 luminosities. 
The distributions of these parameters in the two samples all differ at significance levels exceeding 
99%. These differences are all in agreement with the expectations from hydrodynamical models, 
but for the first time confirmed for a homogeneous and statistically significant sample of galactic 
planetary nebulae. We interpret these differences as evidence for the acceleration of the nebular 
shells during the early evolution of these intrinsically bright planetary nebulae. As is the case 
for planetary nebulae in the Magellanic Clouds, the acceleration of the nebular shells appears to 
be the direct result of the evolution of the central stars. 



Subject headings: 
Galaxy: bulge 

Introduction 



ISM: planetary nebulae (general) — ISM: kinematics and dynamics — stars: evolution- 



work, a variety of theoretical and numerical stud- 
ies of the hydrodynamics have been undertaken to 
understand and predict the kinematic properties 
of planetary nebulae. Our modern view of the 
kinematics of planetary ne bulae was established 
in the early 1990's (e.g.. iKahn fc West I 1 19851: 



Th e interacting stellar winds model o flKwok et al. 
([19781 ) provides a general framework for under- 
standing planetary nebulae. Within this frame- 



1 The observations reported herein were obtained at the 
Observatorio Astronomico Nacional in the Sierra San Pedro 
Martir (OAN-SPM), B. C, Mexico. 



Brcitschwcrdt fc Kahn lll990HKahn fc Breitschwerdt 
19901: iMellema Ill994l ). These studies showed that 
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the AGB envelope is accelerated in two phases. 
The first phase is a result of the shock wave initi- 
ated by the ionization front that sweeps through 
the AGB envelope as the central star's tempera- 
ture increases. Then, as the central star's wind 
energy increases, it drives a pressure-driven central 
bubble that accelerates the AGB envelope through 
ram pressure. More recent numerical work that 
attempts to include more realistic AGB evolution 



confirms these basic results (e.g., IVillaver et al. 
2002HPerinotto et al. |[2004 . Eventually, the cen- 
tral star's wind ceases and the inner part of the 
AGB envelope backfills towards the central star 
while the outer part maintains its momentum - 
driven evolution (e.g-. lGarcia-Segura et al. 1120061 ) . 

These models have been used extensively to 
interpret observational studies of the kinematics 
of many individual planetary nebulae. Unfortu- 
nately, to date, systematic, homogeneous studies 
of planetary nebula populations to compare with 
these theoretical efforts are scarce. iHeap (1993) 
found a correlation between nebular expansion ve- 
locity and the terminal v elocity of the c entra l 
star wind. More recently, iMedina et al. I <|2006h 
claimed a correlation between expansion velocity 
and age indicators such as density and central star 
temperature. However, both of these studies suffer 
from heterogeneity, at least regarding object selec- 
tion. Even regarding planetary nebulae with Wolf- 
Rayet (WR) central stars, the precise role that 
the winds from the central stars play in the kine- 
matics of the nebular shells is not entire ly clear 
(jGesicki et al. 1120061 [Medina et al. 1120061 ). 

The studies of the kinematics of planetary 
nebulae in the Mag ellanic Clouds are perhap s 



the most systematic (jDopita et al. 1 119851 [1988) 



These studies find correlations between expan- 
sion velocity, excitation class, and nebular density, 
suggesting that the properties of the central star 
dominate the nebular evolution. As they point 
out, however, their results are not general since 
they have been found only for a particular popu- 
lation of (bright) planetary nebulae. 

Therefore, there remains a need for a coherent 
population study in a different galactic environ- 
ment. Here, we present a study of the kinemat- 
ics for a large sample of planetary nebulae in the 
bulge of the Milky Way (Bulge). These objects 
were selected in a homogeneous way, with the hope 
of simulating populations of bright extragalactic 



planetary nebulae in bulge-like systems (bulges of 
spiral galaxies, dwarf spheroidals, and elliptical 
galaxies). The observations were all made with the 
same instrument and methodology. Likewise, the 
data reduction and analysis is homogeneous and 
entirely independent of model parameters. Details 
are given in Section 2. We find that the expan- 
sion properties do indeed depend upon the evolu- 
tionary stage of the central star, with planetary 
nebulae hosting hotter central stars having larger 
expansion velocities (Section 3). We also find that 
other parameters that should depend upon age, 
the nebular size, density, and H/3 luminosity, also 
differ significantly between the two groups (Sec- 
tion 4). We compare these properties with theo- 
retical models and find generally good agreement 
(Section 5). We present our conclusions in Section 
6. 

2. Observations and Analysis 

2.1. The Planetary Nebula Sample 

We present the sample of Bulge planetary neb- 
ulae that we observed in Table [TJ There are 
86 objects in total, drawn from existing spec- 



trosco 


pic surveys ( Aller & Keves Ill987c Webster 


1988; 


Cuisinier et al. 1996c iRatag et al. 1 


1997 




Cuisinier et al. 


200C 


: Escudero & Costa 


2001 


Escudero et al. 


2004 


Exter et al. 112004c Gornv et al. 



We selected this sample hoping to simulate 
a sample of bright extragalactic planetary neb- 
ulae in bulge-like systems. Our selection crite- 
ria were refined somewhat over the duration of 
the observations, but eventually converged to (i) 
a position within 10° of the galactic centre, (ii) 
a large observed, reddening-corrected H/3 flux, 
nominally log/(H/3) > — 12.0dcx, (iii) a large 
[O IIIJA5007/H/3 ratio, usually exceeding a value 
of 6, and (iv) the existence of low resolution spec- 
troscopy in which the electron temperature may 
be determined from the [O III]A4363/5007 ratio. 
Usually, the systemic radial velocity exceeds 30 
km/s in an attempt to exclude disk objects to- 
wards the Bulge. We imposed no explicit limit 
on the size of the objects. Nor did we impose an 
upper limit upon the flux (H/3 or 6 cm). In prac- 
tice, however, a flux limit is usually imposed by 
the surveys from which we selected our objects. 

Figure [T] explains the logic of our lower lim- 
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its on the H/3 flux and the [O III]A5007/H/3 ra- 
tio. Basically, bright extragalactic planetary neb- 
ulae in ellipticals, dwarf spheroidals, and the bulge 
of M31 fall within these limits (see Figure cap- 
tion). Restricting the sample to objects with large 
[O III]AA5007/H/3 ratios will exclude planetary 
nebulae very early in their evolution when their 
central stars are coolest. Our limit on the H/3 lu- 
minosity does not restrict the planetary nebulae 
to contain central stars that are on the horizon- 
tal portion of their evolutionary track, but also 
includes the initial fading following the extinction 
of nuclear reactions. 

2.2. Observations and Data Reductions 

We acquired our observations during eight ob- 
serving runs spanning the period from 2003 June 
to 2007 August at the Observatorio Astronomico 
Nacional in the Sierra San Pedro Martir, Baja Cal- 
ifornia, Mexico (OAN-SPM). More details of the 
observations will be provided in a forthcoming pa- 
per. 

High resolution spectra were obtained with 
the Manchester echelle spec trometer (MES-SPM; 
iMeaburn et al. Hl984l l2003fh The MES-SPM is a 
long slit, echelle spectrometer that has no cross- 
dispersion. Instead, narrow-band filters were used 
to isolate orders 87 and 114 containing the Ha 
and [O III]A5007 emission lines, respectively. A 
150 ^m wide slit was used for the observations, 
resulting in a slit l'/9 wide on the sky. Coupled 
with a SITe 1024 x 1024 CCD with 24 ^m pixels 
binned 2x2, the resulting spectral resolutions 
were approximately 0.077A/pix and O.lOOA/pix 
at [O III]A5007 and Ha, respectively (equivalent 
to llkm/s for 2.6pix FWHM). The spectra were 
calibrated using exposures of a ThAr lamp taken 
immediately before or after every object exposure. 
The internal precision of the arc lamp calibrations 
is better than ±1.0km/s. 

Usually, only one deep spectrum was obtained 
in each of the [O III]A5007 and Ha filters. The 
exposure times varied depending upon the bright- 
ness of the object. At most, the exposure times 
were of 30 minutes duration. If this was expected 
to produce saturation, shorter exposure times were 
used. The exposure time for the Ha spectrum was 
chosen to achieve a signal level similar to that in 
the deep [O III]A5007 spectrum, based upon the 
intensities of these lines observed in the low reso- 



lution spectra, but was also limited to a maximum 
of 30 minutes duration. 

All of the planetary nebulae observed are re- 
solved (see Table [1]). In all cases, the slit was cen- 
tered on the object. Usually, the slit was oriented 
in the north-south direction. 

The spectra were reduced using the twodspec 
and specred packages of the Image Reduction and 
Analysis FacilitjH (IRAF). The pro cedure for data 



reduction followed that outlined in Massev et al 



( 1992L Appendix B) for long slit spectroscopy. The 
object spectra were edited of cosmic rays. Then, 
a nightly mean bias image was subtracted from 
each object spectrum. Next, the arc lamp spec- 
tra were used to map positions of constant wave- 
length. These maps were then used to rectify the 
object spectra so that lines of constant wavelength 
fell exactly along the columns, a process that si- 
multaneously applies a wavelength calibration. 
Finally, wavelength-calibrated, one-dimensional 
spectra were extracted for each object. No flux 
calibration was performed. 

2.3. Kinematic parameters 

We shall present all of the one-dimensional line 
profiles in a subsequent publication. For that rea- 
son, we do not present the line profiles here. In- 
stead, we turn to an explanation of the analysis of 
the kinematics. 

The one-dimensional spectra were analyzed us- 
ing a l ocally-implemented software package (IN- 
TENS; iMcCall et al"lll985l ) to determine the ra- 
dial velocity, flux, and profile width (FWHM; 
full width at half maximum intensity) as well as 
the uncertainties in these parameters. This soft- 
ware fits the emission line profile with a sampled 
gaussian function and models the continuum as a 
straight line. Thus, this analysis assumes that the 
lines have a gaussian shape and that they are su- 
perposed on a flat continuum. In the case of the 
Ha line, the He IIA6560 line may also be present. 
In this fit is made simultaneously to both 

lines and the continuum. 

While the assumption of a gaussian line shape 
is reasonable for bright extragalactic planetary 



2 IRAF is distributed by the National Optical Astronomical 
Observatories, which is operated by the Associated Uni- 
versities for Research in Astronomy, Inc., under contract 
to the National Science Foundation. 
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nebulae, because they are spatially unresolved, 
it might seem odd for spatially-resolved objects. 
Nonetheless, for the large majority of the objects, 
the one-dimensional line profiles are usually not 
double-peaked and the deviations from a gaussian 
shape usually represent less than 1 0% of the to- 



tal fl ux (see Fig. [2] for an example; iRicher et al. 
2009h . 



Table [T] presents the observed line widths 
(FWHM) and their uncertainties for each object 
in both Ha and [O III]A5007. The uncertainties in 
the line widths are the formal uncertainties (one 
sigma) from fitting a sampled gaussian function 
to the line profile. In order to derive the intrin- 
sic line widths, the observed line widths must be 
corrected for several effects, all of which are as- 
sumed to contribute to the observed line width 
in quadrature. The effects that broaden the true, 
intrinsic profile arc instrumental (<7j nst ), thermal 
(<7th), and fine structure (cr/ s ) broadening, 



obs 



The first term, ol. 



'th 



(1) 



is the true, intrinsic line 
width resulting from the kinematics of the plan- 
etary nebula. The instrumental profile has a 
measured FWHM of 2.5-2.7 pixels. We adopt a 
FWHM of 2.6 pixels for all objects, which amounts 
to about 11 km/s (FWHM). We compute th e ther- 



mal broadening from the usual formula (jLang 
1980l eq. 2-243), adopting rest wavelengths of 



6562.83A and 5006.85A for Ha and [O III]A5007 
respectively, and assuming no turbulent veloc- 
ity. The thermal broadening (FWHM) amounts to 
0.47A (21.4km/s) and 0.089A (5.3km/s) for Ha 
and [O III]A5007, respectively. The fine structure 
broadening was taken to be <j f s = 3. 199 km/s for 
Ha and zero for [O III]A5007 (|Garci'a-Diaz et al. 
2008h . 



The interpretation of the resulting line width. 

Ay, 



AV = 2.3556cr t „, 



(2) 



is not necessarily simple. In the case of an expand- 
ing, spherical shell that is spatially-unresolved, the 
line width corresponds to the expansion velocity. 
If the velocity field or matter distribution is more 
complicated or if the object is resolved spatially, 
then the line width is a luminosity-weighted veloc- 



ity width for the mass projected within the spec- 
trograph slit. 

The velocity width that we measure should typ- 
ically exceed the luminosity-weighted line width 
for the entire object. The spectrograph slit was 
centered on each object and all of the objects are 
wider than the slit (Table [T]). Therefore, matter 
near the edges of the objects is excluded and this 
matter is likely to have the lowest projected veloc- 
ity along the line of sight. Since the matter at the 
edges of the objects would contribute low system- 
atic velocities for each object, the line profile we 
measure for the matter included within the slit will 
be slightly larger than the true mass-weighted line 
width. These argu ments are supporte d by th e re- 
sults presented by iGesicki fc Ziilstra" (2000) and 



Rozas et al. I (|2007h . Their simulations of thin, 
expanding, spherical shells indicate that the line 
widths we measure may over-estimate the inte- 
grated line widths for the entire objects by up to 
approximately 15%, but that the exact amount 
will depend upon the fraction of the object cov- 
ered by the slit. 

Fortunately, the discussion that follows does 
not depend upon any interpretation of the line 
width. Since it is clear, however, that this line 
width will be similar to twice the expansion ve- 
locity, Table [T] presents half of the line width in 
velocity units for each object, i.e., 



AVb.5 = 0.5AV = 1.1778cr tr 



(3) 



3. Nebular Kinematics versus the Evolu- 
tionary State of the Central Star 

For about half of the objects in our sample, the 
He II A6560 line is present in the Ha spectra. Fig. 
[2] presents M 2-16 as an example. In practice, 
we can detect He II A6560 provided its intensity 
is at least 0.4% that of Ha and a signal-to-noise 
of at least 50 is achieved for the Ha line. For 
comparison with theoretical models, this limit is 
predicted to occur for an effective temperature of 
about 9 x 10 4 K acc ording to the hydrodynam ical 
models presented in ISchonberner et al. I (|20071 ). if 
we use the the atomic dat a for hydrogen and he- 
lium in lOsterbrockl (|l989h . Table [2 indicates the 
objects for which the He II A6560 line was observed 
in our Ha spectra. 

Obviously, the objects in which the He II A6560 
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line appears are those with the hotter central stars. 
Therefore, we may separate our sample into two 
samples with cooler and hotter central stars ac- 
cording to whether the He IIA6560 line is ob- 
served. Given our selection criteria, this separa- 
tion into samples with cooler and hotter central 
stars is equivalent to a separation into samples in 
different stages of evolution. The planetary neb- 
ulae with the hotter central stars are in a more 
advanced stage of evolution than those with the 
cooler central stars. 

Fig. [3] presents histograms of the line width 
(AV0.5) distributions for the two samples of plan- 
etary nebulae. Clearly, the distributions appear 
different. The sample of planetary nebulae with 
the cooler central stars have predominantly low 
line widths, while the sample with the hotter cen- 
tral stars have a much broader distribution of line 
widths that extends to substantially higher values. 

A variety of statistical tests may be applied 
to these distributions to determine whether it is 
probable that they may result from a single par- 
ent distribution. Since we do not know a priori 
the form of the distribution that the line widths 
should have, non-parametric statistical tests are 
most appropriate. Perhaps, the best-known of 
these tests is the Kolmogorov-Smirnov two-sample 
test to the cumulative distribution functions of 
the line widths. However, a more powerful non- 
parametr ic test is the Wilcoxon-M ann- Whitney li- 
test fe.g.. lWall fc Jenkins Il2003l §5.4.3). This is a 
rank test in which the line widths from the two 
samples are combined into a single sample and 
ranked. Then, the rankings of the line widths for 
the two samples are compared with a uniform dis- 
tribution of rankings. 

The U-test indicates that the probability of 
drawing the two samples of line widths from the 
same parent distribution is only 1.3 x 10~ 7 . The 
nebular shells surrounding hotter central stars are 
expanding systematically more rapidly than their 
counterparts surrounding cooler central stars. 
Since all of these hotter central stars were nec- 
essarily cooler in the past, the simplest conclusion 
is that the nebular shells in bright planetary neb- 
ulae are accelerated as the star gets hotter. 



4. Other Properties versus the Evolution- 
ary State of the Central Star 

Considering that the nebular shells of bright 
planetary nebulae are accelerated as the central 
star gets hotter, one might expect other proper- 
ties of the two samples of planetary nebulae to 
differ. Essentially, characteristics that vary sys- 
tematically with the age or evolutionary stage of 
the planetary nebula should differ between the two 
samples. 

The most obvious characteristic to consider is 
the nebular size. The cooler central stars should 
be found in younger planetary nebulae that have 
had less time to expand, and should be system- 
atically smaller than those surrounding the hot- 
ter central stars. There are no systematic mea- 
surements of the diameters of the planetary neb- 
ulae in our sample, so we measured nebular di- 
ameters from our deep Ha spectra. To determine 
the diameters, we collapsed the spectra along the 
wavelength axis to produce one-dimensional spa- 
tial profiles. We then measured the nebular sizes 
at 50% and 10% of the peak intensity. These di- 
ameters are found in Table [TJ 

The distributions of diameters at 10% of the 
peak flux are shown in Fig. |U Again, it is clear 
that the sample with the hotter central stars has 
a distribution of diameters that tends to larger 
sizes than the sample of planetary nebulae with 
the cooler central stars. Applying the U-test, we 
obtain a probability of only 7 x 10~ 4 that the 
two samples arise from a single parent population. 
Hence, the nebulae surrounding the hotter central 
stars are larger than those surrounding the cooler, 
less evolved, central stars. 

Conceivably, the larger sizes of the nebulae 
around hotter central stars might artificially bias 
their line widths to larger values because the fixed 
spectrograph slit covers a smal ler fraction of these 



larger neb ulae. As shown by Gesicki &: Ziilstra 



( 2000f ) and Rozas et al. ( 2007 ). this might cause 
the line widths for larger nebulae to be over- 
estimated. To test whether this is the case, we 
estimated how much the line widths could be over- 
estimate d based upon the re s ults p resented in Ta- 



ble 1 of Gesicki & Ziilstra (2000). We then di 



vided the line widths of all objects in Table Q] 
(AV0.5) by this correction and re-calculated the 
probability that the two samples arise from a sin- 
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gle parent population. The probability increases, 
but to only 2.5 xlCT 7 (considering the [O III]A5007 
line widths), so the conclusion that the two distri- 
butions are statistically very different is not af- 
fected. 

The nebular density is also expected to evolve, 
decreasing as the nebular volume increases. Table 
Q] lists the electron densities determined from the 
[S II]AA6716,6731 ratio as measured in low resolu- 
tion spectroscopy from the literature. When there 
was more than one measurement, the measure- 
ments were averaged. The distributions of neb- 
ular density in the two samples are shown in Fig. 
[5l In this case, the difference is more subtle, but 
the planetary nebulae with the hotter central stars 
tend to have lower densities. Applying the U-test, 
we find a probability of 8 x 10 -3 that the two 
samples arise from the same parent population. 
Therefore, the planetary nebulae with hotter cen- 
tral stars have systematically lower electron den- 
sities. 

Finally, we compare the H/3 luminosities of the 
two samples. We expect all of the planetary nebu- 
lae without He IIA6560 emission to contain central 
stars on the horizontal portion of their evolution- 
ary tracks. It is unlikely that any of these objects 
contain central stars that are fading towards the 
white dwarf regime since virtually all of the cen- 
tral stars should be sufficiently massive to achieve 
temperatures in excess of 10 5 K before fading, by 
which time emission from He 2+ ions should be ev- 
ident. On the other hand, the sample of planetary 
nebulae with hotter central stars should contain 
objects whose central stars are fading towards the 
white dwarf regime. Therefore, we would expect 
that the H/3 luminosities of the two samples could 
differ, with the sample with the cooler central stars 
having higher luminosities. 

The H/3 luminosities for the Bulge planetary 
nebulae in Tableland Fig. [T]are not very accurate 
for a variety of reasons. First, the observed H/3 
fluxes are usually those measured through a spec- 
trograph slit and not from photometry. Likely, 
the H/3 fluxes are underestimated for many Bulge 
planetary nebulae. Second, these fluxes must be 
corrected for reddening. The reddening may in- 
clude instrumental or observational effects, such 
as differential atmospheric refraction, and so may 
over-estimate the true reddening. These redden- 
ings have been used to correct the H/3 fluxes for 



extinction. Third, the individual distances for the 
planetary nebulae are unknown. To compute the 
H/3 luminosities in Table [TJ we adopted a common 
distance of 7.5 kpc. In all likelihood, the mean 
distance for the sample will be less than this and 
it is probably this assumption that leads to larger 
H/3 luminosities than are observed in extragalactic 
planetary nebulae. 

Fig. [5] presents histograms of the H/3 luminos- 
ity distributions for the two samples, divided ac- 
cording to the presence or absence of He IIA6560 
emission. The H/3 luminosities are taken from Ta- 
ble Q] Clearly, the distribution for the sample with 
the cooler central stars extends to higher H/3 lumi- 
nosities. The mean luminosities for the two sam- 
ples differ by 0.36 dex, or 0.90 mag. The U-test 
indicates that the probability that two distribu- 
tions arise from the same parent population is only 
2.2 x 10~ 4 . Therefore, the planetary nebulae with 
the cooler central stars are brighter, on average, 
than those with the hotter central stars. 

Summarizing our results, we find significant dif- 
ferences between our samples of planetary nebulae 
with cooler and hotter central stars. The plane- 
tary nebulae with cooler central stars have smaller 
line widths, smaller sizes, higher densities, and 
higher H/3 luminosities. The differences are all sig- 
nificant at confidence levels exceeding 99% (Table 
©■ 

5. Discussion 

At present, the theory of the evolution of plan- 
etary nebulae outlined earlier is unable to predict 
definitive expansion velocities for planetary nebu- 
lae or the AGB stars from which they arise. Obser- 
vations of AGB stars find that their winds usually 
have an expansion velocit y below 15km/s and of- 
ten below l Okm/s (e.g., iRamstedt et al. 1 120061: 
~ 199lh . Observations of planetary nebulae 



Lewis 



typical ly find expansion velociti es of 15-30 km/s 



Gesicki fc Ziilstral l2000t iMedina et at 



2006). At least qualitatively, theory may explain 
the difference: The ionization front drives a shock 
wave through the AGB envelope, raising the ex- 
pansion velocit y by ~ 5 — 6 km/s over that of the 



AGB envelope ([Chevalier Ill997t iPerinotto et al. 
20041 ). Then, the pressure-driven bubble expands 



into the the AGB envelope, sweeping up a dense, 
inner rim of matter that may increase the ex- 
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pansion velocity significantly, particularly if this 
rim is able to overtake the density enhancement 
produced by the shock wave. 

The median line width in [O III]A5007 (Ha) for 
our sample of Bulge planetary nebulae with cooler 
central stars is 17.2 km/s (15.4 km/s) and the 
standard deviation of the distribution is 5.1 km/s 
(4.0 km/s). Since we expect our line widths to 
be slightly larger than the expansion velocity, the 
median line width for the sample with the cooler 
central stars is approximately that expected if the 
AGB envelopes had expansion velocities of about 
10 km/s, as observations suggest commonly occurs 



(Mj 5 2M^: e.g.jMellema 



199 1 IVillaver" et al. 



( Ramstedt et al. 20061 ). It is therefore natural 



to associate these planetary nebulae with models 
during the stage in which the ionization front is 
driving a shock wave through the AGB envelopes. 

The median line width in [O III]A5007 (Ha) 
for the sample of Bulge planetary nebulae with 
the hotter central stars is substantially higher, 
23.7 km/s (24.3 km/s), as is the dispersion in ve- 
locities, 7.6 km/s (7.1 km/s). In the models, by 
the time that the central star has a tempera- 
ture of 9 x 10 4 K, the stellar wind-powered bubble 
is expanding into the ionized AGB envelope and 
raising the observed expansion velocities. Thus, 
our sample of planetary nebulae with hotter cen- 
tral stars is naturally associated with this sec- 



ond stage of envelope acceler ation (|Villaver et al. 
20021; IPerinotto et al. |[200i ). 



It is more difficult to compare the differences 
in size that we find for our two samples with 
existing theoretical models. In part, this has 
to do with how we have determined our nebu- 
lar sizes. The main problem is the mismatch in 
spatial resolution between our observations and 
one-dimensional theoretical profiles. Our obser- 
vations generally include a substantial fraction of 
the object, while one-dimensional theoretical pro- 
files have infinite spatial resolution. In spite of 
this difficulty, the available theoretical profiles in- 
dicate that the outermost dense/b right regions in- 
crease in size as time progresses (iMelle ma 111994; 



Villaver et a"l~ll2002t IPerinotto et al. 112004 ). even 



though the size of the brightest str ucture may not 



incre ase monotonically with time (|Villaver et al 
20021) . 

Theoretical models indicate that the entire 
AGB envelopes should become ionized relatively 
early for the stellar masses of relevance here 



2002t IPerinotto et al. 1 120041 ). They are also ex- 
pected to remain entirely ionized. Once the en- 
velope is entirely ionized, its H/3 luminosity will 
decrease with tim e as a result of dilution (e.g., 
Osterbrock I ll989l) . As a result, that we find lower 



H/3 luminosities for the sample of objects that is 
more evolved is entirely compatible with theoreti- 
cal expectations. 

Likewise, theoretical models are compatible 
with our findings concerning the electron density. 
Theoretical models indicate that the densities are 
higher while the expansion velocities are domi- 
nated by the ionization front than later on when 
the pressure-driven bubble dominates. However, 
this agreement is bitter-sweet. We use densities 
derived from [S II] lines. In principle, these lines 
from the recombination zone should not exist in 
the objects with the hotter central stars, since 
these objects are predicted to have completely 
ionized their AGB envelopes. In practice, this 
emission no doubt arises in density enhancements 
that are absent from the models, and so the [S II] 
densities are likely an upper limit to the typical 
densities in these objects. Nonetheless, the mod- 
els predict a general dilution of the entire AGB 
envelope and this dilution is probably what gives 
rise to the decay in [S II] densities that we observe. 

It is not unusual that the line width distribu- 
tions overlap for our samples of planetary nebu- 
lae with cooler and hotter central stars. As the 
theoretical models make very clear, the expansion 
velocity depends upon the AGB envelope expan- 
sion velocity and the mass loss rate at the end of 
the AGB evolution. Undoubtedly, these properties 
varied among the stellar progenitors of the plane- 
tary nebulae in both samples. Furthermore, even 
if they did not, the planetary nebulae in both sam- 
ples are found in a range of evolutionary states, 
particularly for the sample with the hotter central 
stars. Consequently, the line widths we observe 
will have varying contributions from the acceler- 
ations due to the shock front and the pressure- 
driven bubble. Unless the objects are all spheri- 
cal (unlikely), the different orientations will also 
broaden the line width distributions. 

It is also not unusual that the distributions of 
line widths differ more than do those for the neb- 
ular diameters, densities, and luminosities. The 
structure that appears brightest varies with time 
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(jVillaver et al. 1120021 ; [Perinotto et aT~l 12004 ). de- 



pending upon the development of the density en- 
hancements due to the ionization front or the 
swept-up inner rim. It is therefore not surpris- 
ing that there is less difference between the dis- 
tributions of densities and diameters than there 
is between the distributions of line widths. Like- 
wise, since the H/3 luminosity will depend upon 
how much mass is sw ept up by the AGB envelope 
(|Villaver et al. Il2002h . and that this may vary sig- 
nificantly from one object to the next depending 
upon location, the H/3 luminosities may be some- 
what scattered. 



Our findings confir m those of lHeap I (| 19931) and 
Medina et al. I (|2006h . who found correlations be- 
tween the nebular kinematics and the evolutionary 
state of the central star. Our findings, however, 
are based upon a sample of planetary nebulae that 
is larger, more homogeneous, and whose analysis 
is entirely independent of models. It is perhaps 
not surprising that our results demonstrate much 
more clearly the intimate relationship between the 
central star's evolution and the kinematics of the 
nebular shell. 

Perhaps, the m ost relevant st u dies for com par- 
ison are those of Dopita et al. ( 19851 19881 ) for 



planetary nebulae in the Magellanic Clouds. They 
found that the nebular expansion velocity depends 
upon the excitation class and H/3 luminosity of the 
nebula. Both the excitation class and the H/3 lu- 
minosity are functions of the evolutionary stage 
of both the central star and nebular shell. In 
other words, the kinematics of planetary nebulae 
in the Magellanic Clouds require a coordination 
between the evolutionary states of the central star 
and the nebular shell. Our results clearly demand 
the same sort of coordination between the evolu- 
tion of the central star and the surrounding neb- 
ular shell. 

The exact relationship between the central star 
and the nebular kinematics need not be conserved 
between the Magellanic Clouds and the Bulge. In- 
deed, it is difficult to compare the results directly 
as many complications may arise. Ideally, the neb- 
ular kinematics should be related to the properties 
of the central star in both environments (temper- 
ature, wind velocity, mass loss rate). However, for 
most of the objects in both the Bulge and the Mag- 
ellanic Clouds, observations of the central star's 
wind are unavailable. Excitation classes or line ra- 



tios provide indirect temperature diagnostics, but 
these will usually depend upon the chemical com- 
position of the nebula and are also susceptible to 
systematic biases if there are differences in the dis- 
tributions of the central star masses or AGB wind 
properties. 

(l2006f ). 



Our findings, like those of lGesicki et al 



do not support the claim that the nebular shells 
surrounding WR central stars have larger expan- 
sion velocities ( Medina et al. II2OO6 ). We made no 
special attempt to include planetary nebulae with 
WR central stars in our sample, so it includes only 
eight examples (see Table [lj. Of these, five of the 
central stars are cool and three are hot enough to 
produce He II emission. With such small numbers, 
the statistics are clearly not conclusive. While the 
cooler WR central stars are found in nebulae with 
a higher mean velocity than the rest of the plan- 
etary nebulae with cool central stars, the hotter 
WR central stars are found in nebulae with the 
same line widths as the rest of the sample of hot- 
ter central stars. Perhaps, WR central stars are 
able to initiate the wind-driven phase of nebular 
acceleration earlier, but they may not necessarily 
afford greater acceleration in the long run. The 
planetary nebulae with normal central stars in 
our sample span a larger range in line width than 
do the expansion velocities obser ved in nebulae 
surrou nding WR central stars by iMedina et al. 
(2006). Resolving this issue will require a more 
careful, dedicated study. 

Therefore, our observations of bright planetary 
nebulae in the Milky Way bulge confirm the kine- 
matics predicted by theoretical models. It is no- 
table that observations not only recover the se- 
quence of nebular acceleration that is predicted, 
but also the coincidence of these stages with par- 
ticular properties of the central stars. The details 
of the models differ considerably in parameters 
that are difficult to both observe and model, such 
as the structures within the AGB envelopes or the 
transition time between the AGB and planetary 
nebula stages. 

6. Conclusions 

We have obtained kinematic data for a large 
sample of planetary nebulae in the Milky Way 
bulge, selected with the goal of simulating sam- 
ples of bright extragalactic planetary nebulae in 



S 



bulge-like environments. For most of the sample, 
our criteria also included a reddening-corrected 
H/3 flux exceeding logF(H/3) > — 12.0dex and 
[O III]A5007/H/3 > 6, in addition to the usual 
criterion of projected proximity to the Galac- 
tic center. We measure line widths for Ha and 
[0 III]A5007. Our sample is the largest that has 
been selected and observed in a homogeneous way. 
We have also analyzed it in a completely model- 
independent manner. 

For half of the sample, the He IIA6560 line is ob- 
served in the Ha spectra. Since this line appears 
only for hotter central stars, it allows us to divide 
our sample according to the evolutionary stage of 
the central star, with the more evolved objects 
found in the subsample containing the hotter cen- 
tral stars. The kinematics of the two samples are 
significantly different, with the sample containing 
the hotter central stars expanding faster. We also 
compare the diameters, electron densities (from 
[S II] AA6716, 6731), and H/3 luminosities for the 
two subsamples. In all cases, there are statisti- 
cally significant differences. The planetary nebu- 
lae hosting the hotter central stars are larger, less 
dense, and less luminous than their counterparts 
with cooler central stars. All of these differences 
exceed a statistical significance of 99%. Also, all 
of these differences are compatible with the results 



of hydrodynamical models (e.g., Me 


ema 1994 


Villaver et al. 2002: Perinotto et al. 


2004). 



Our primary conclusion is that we have clearly 
observed the acceleration of the nebular shells in 
planetary nebulae in the bulge of the Milky Way 
and that this occurs during the early evolution 
of their central stars. Our findings, based upon 
a large, homogeneous sample, constitute the first 
unequivocal evidence that the nebular kinematics 
depend upon the evolutionary state of the cen- 
tral star for planetary nebulae in the Milky Way 
and are far clearer and more convincing than any 
previous results. Hence, there is now very clear 
evidence that this dependence occurs in at least 
two environments: the Magellanic Clouds and the 
bulge of the Milky Way. Given the differences in 
the stellar populations in these two environments, 
it is likely that the dependence of the nebular kine- 
matics on the evolutionary state of central star 
is more general and that it occurs in all environ- 
ments. 
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Table 1 

Bulge Planetary Nebula Sample 



Object PNG Run FWHM(Ha) 1 AV .5(Ha) FWHM(5007) a AV . 5 (5007) diameter (in Ha) b N e " log F(H/3) 6560? WR? d 











( A ) 




(Km / s ) 


( A ) 




^Km/s ) 




OU/0 l m ax 




(cm ) 


(crg/s) 






Bl 3-13 


000.9-02.0 


2006 


Jun 


0.9427 ± 


.0028 


16.05 ± 0.06 


0.6063 ± 


.0079 


16.93 ± 0. 


.24 


2.4 


5.5 


1765 


34.19 


no 


wcls 


Cn 1-5 


002.2-09.4 


2004 


Jun 


1.0384 ± 


.0026 


18.87 ± 0.06 


1.251 ± 0. 


.041 


36.9 ± 1. 


2 


4.4 


8.6 




35.21 


no 


WR 


Cn 2-1 


356.2-04.4 


2004 


Jun 


0.9255 ± 


.0039 


15.51 ± 0.09 


0.6826 ± 


.0074 


19.36 ± 


.22 


2.7 


4.8 


6967 


34.97 


no 


no 


H 1-1 


343.4+11.9 


2004 


Jun 


1.617 ± 


.019 


34.04 ± 0.43 


1.238 ± 


.043 


36.5 ± 1 


.3 


2.5 


4.4 




33.88 


yes 


no 


H 1-11 


002.6+08.2 


2006 


Jun 


0.9219 ± 


.0057 


15.40 ± 0.13 


0.625 ± 


Oil 


17.53 ± 0. 


.32 


4.4 


7.3 


10803 


34.43 


yes 


wcls 


H 1-14 


001.7+05.7 


2005 


Jul 


1.803 ± 


.036 


38.62 ± 0.81 


1.313 ± 


.068 


38.8 ± 2 





4.3 


7.4 


985 


34.30 


yes 


no 


H 1-16 


000.1+04.3 


2005 


May 


1.026 ± 


.003 


18.51 ± 0.07 


0.743 ± 


.012 


21.25 ± 


36 


2.8 


4.7 


7035 


35.11 


yes 


no 


H 1-17 


358.3+03.0 


2005 


Jul 


0.9635 ± 


.0035 


16.67 ± 0.08 


0.6150 ± 


.0072 


17.20 ± 


22 


2.1 


4.7 


16190 


35.29 


no 


no 


H 1-18 


357.6+02.6 


2004 


Jun 


0.8105 ± 


.0026 


11.67 ± 0.06 


0.5062 ± 0. 


.0039 


13.66 ± 


.12 


2.6 


4.8 


8627 


34.94 


yes 


no 


H 1-20 


358.9+03.2 


2003 


Jun 


0.9361 ± 


.0025 


15.85 ± 0.06 


0.6905 ± 


.0042 


19.53 ± 


13 


2.4 


4.9 


3811 


34.85 


no 


no 


H 1-23 


357.6+01.7 


2005 


May 


1.0266 ± 


0080 


18.54 ± 0.18 


0.754 ± 


016 


21.62 ± 


.49 


3.2 


5.9 


3840 


35.05 


no 


no 


H 1-27 


005.0+04.4 


2003 


Jun 


1.0238 ± 


.0029 


18.48 ± 0.07 


0.7004 ± 


.0075 


19.94 ± 


22 


1.7 


3.8 


18385 


34.83 


no 


no 


H 1-30 


352.0-04.6 


2006 


Jun 


0.9616 ± 


.0044 


16.62 ± 0.10 


0.6000 ± 


.0052 


16.73 ± 0. 


.16 


3.2 


4.4 


6495 


34.20 


yes 


no 


H 1-31 


355.1-02.9 


2005 


May 


0.9932 ± 


.0037 


17.56 ± 0.08 


0.6646 ± 


.0038 


18.79 ± 0. 


.11 


2.8 


5.1 


12957 


34.84 


yes 


no 


H 1-32 


355.6-02.7 


2005 


May 


0.7685 ± 


.0026 


10.08 ± 0.06 


0.4385 ± 


.0033 


11.37 ± 0. 


.10 


3.1 


5.2 


7897 


34.04 


no 


no 


H 1-33 


355.7-03.0 


2004 


Jun 


0.7904 ± 


.0018 


10.93 ± 0.04 


0.5218 ± 


.0022 


14.17 ± 


.07 


3.9 


6.5 


3215 


34.72 


no 


no 


H 1-40 


359.7-02.6 


2005 


May 


1.014 ± 0. 


.014 


18.16 ± 0.32 


0.643 ± 


.014 


18.09 ± 


.42 


3.1 


7.0 


11745 


35.42 


no 


no 


H 1-41 


356.7-04.8 


2006 


Jun 


1.1607 ± 


.0091 


22.30 ± 0.21 


0.832 ± 


.014 


24.05 ± 


.43 


5.7 


8.8 


2525 


34.35 


yes 


wcls 


H 1-42 


357.2-04.5 


2006 


Jul 


0.9351 ± 


.0059 


15.81 ± 0.13 


0.5362 ± 


.0080 


14.65 ± 


.24 


3.0 


5.2 


2270 


34.92 


no 


wcls 


H 1-45 


002.0-02.0 


2005 


Jul 


1.590 ± 


.011 


33.37 ± 0.26 


1.016 ± 


.011 


29.70 ± 


.34 


2.7 


4.4 




35.03 


yes 


no 


H 1-50 


358.7-05.2 


2004 


Jun 


1.0263 ± 


.0018 


18.53 ± 0.04 


0.7291 ± 


.0038 


20.82 ± 


.11 


2.2 


4.5 


6920 


34.95 


yes 


no 


H 1-54 


002.1-04.2 


2007 


Aug 


0.9029 ± 


.0045 


14.81 ± 0.10 


0.5304 ± 


.0028 


14.50 ± 


08 


2.8 


5.2 


11594 


34.78 


no 


no 


H 1-56 


001.7-04.6 


2007 


Aug 


0.8261 ± 0. 


.0013 


12.25 ± 0.03 


0.5650 ± 


.0048 


15.67 ± 


.14 


3.2 


6.1 


1164 


34.52 


no 


wcls 


H 1-59 


003.8-04.3 


2005 


May 


1.237 ± 0. 


.015 


24.33 ± 0.34 


1.115 ± 0. 


.033 


32.76 ± 


.99 


3.6 


6.8 


1100 


34.26 


yes 


no 


H 1-60 


004.2-04.3 


2005 


May 


0.9560 ± 


.0076 


16.45 ± 0.17 


0.672 ± 0. 


.013 


19.04 ± 0. 


.38 


3.2 


7.0 




34.26 


no 


wcls 


H 1-67 


009.8-04.6 


2005 


Jul 


1.474 ± 


.030 


30.46 ± 0.68 










5.2 


9.5 


1455 


34.36 


yes 


WR 


H 2-10 


358.2+03.5 


2004 


Jun 


1.1028 ± 


.0093 


20.70 ± 0.21 


0.7907 ± 0. 


.0038 


22.75 ± 0. 


.11 


1.9 


4.1 


8140 


34.61 


no 


no 


H 2-11 


000.7+04.7 


2005 


Jul 


0.7697 ± 


.0058 


10.13 ± 0.13 


0.4684 ± 


.0055 


12.39 ± 


16 


2.6 


5.2 


15550 


35.03 


no 


wcls 


H 2-18 


006.3+04.4 


2004 


Jun 


1.473 ± 0. 


.020 


30.44 ± 0.47 


1.118 ± 0. 


.028 


32.84 ± 


.85 


2.9 


6.6 




34.14 


yes 


no 


Hb 8 


003.8-17.1 


2004 


Jun 


0.9662 ± 


.0022 


16.75 ± 0.05 


0.6465 ± 0. 


.0037 


18.21 ± 


.11 


1.8 


4.1 


14420 


34.05 


no 


no 


He 2-250 


000.7+03.2 


2003 


Jun 


1.185 ± 


.021 


22.96 ± 0.48 


0.797 ± 


.013 


22.98 ± 


40 


3.6 


7.9 


2430 


34.62 


yes 


no 


Hf 2-1 


355.4-04.0 


2005 


May 


1.774 ± 


.059 


37.9 ± 1.3 


1.285 ± 


.069 


37.9 ± 2. 


.1 


9.1 


19.7 




34.17 


yes 


no 


K 5-1 


000.4+04.4 


2006 


Jul 


0.8709 ± 0. 


.0027 


13.76 + 0.06 


0.6191 ± 


.0078 


17.34 ± 


.23 


3.6 


5.7 


1400 


34.61 


no 


wels 


K 5-3 


002.6+05.5 


2006 


Jul 


1.368 ± 


.041 


27.77 ± 0.94 


0.538 ± 


.020 


14.70 ± 0. 


.59 


2.6 


9.5 


1300 


34.01 


yes 


WR 


K 5-4 


351.9-01.9 


2006 


Jul 


0.8187 ± 


.0013 


11.97 ± 0.03 


0.4946 ± 


.0016 


13.28 ± 


.05 


2.2 


3.4 


7600 


34.99 


no 


no 


K 5-5 


001.5+03.6 


2006 


Jul 


0.9265 ± 0. 


.0047 


15.54 ± 0.11 


0.612 ± 


.011 


17.12 ± 0. 


.31 


2.5 


4.1 


10900 


34.85 


no 


no 


K 5-6 


003.6+04.9 


2006 


Jul 


1.302 ± 


.018 


26.04 ± 0.42 


0.866 ± 


.027 


25.09 ± 


.82 


3.6 


8.6 


1000 


33.85 


yes 


no 


K 5-7 


003.1+04.1 


2006 


Jul 


1.619 ± 


023 


32.49 ± 0.52 


1.323 ± 


035 


39.1 ± 1 


.0 


8.9 


13.2 


300 


33.62 


yes 


no 


K 5-9 


355.54-1.4 


2006 


Jul 


1.289 ± 


.035 


25.72 ± 0.81 


0.964 ± 


.043 


28.1 ± 1. 


3 


4.0 


7.7 


3500 


34.60 


yes 


no 


K 5-11 


002.3+02.2 


2006 


Jul 


0.9577 ± 


.0054 


16.51 ± 0.12 


0.6465 ± 


.0087 


18.22 ± 


.26 


5.6 


10.8 


604 


34.31 


no 


no 


K 5-12 


353.5-03.3 


2006 


Jul 


1.548 ± 0. 


.017 


32.33 ± 0.38 


1.208 ± 


.031 


35.59 ± 0. 


.93 


4.1 


7.3 


4480 


34.07 


yes 


no 


K 5-14 


003.9+02.6 


2007 


Aug 


1.0314 ± 


.0087 


18.69 ± 0.20 


0.5986 ± 0. 


.0052 


16.50 + 


.16 


2.0 


3.6 


8000 


34.16 


yes 


no 


K 5-17 


004.3+02.1 


2007 


Aug 


1.435 ± 


.027 


29.48 ± 0.62 


1.043 ± 


.037 


30.5 ± 1 


.1 


2.4 


5.1 


9960 


34.41 


yes 


no 


K 5-19 


005.1+02.0 


2007 


Aug 


1.409 ± 


016 


28.81 ± 0.36 


1.191 ± 


030 


35.08 ± 


.91 


3.0 


5.9 


1170 


34.15 


yes 


no 



Table 1 — Continued 



Object 


PN G 


Run 


FWHM(H 




AVo. 5 (Ha) 


FWHM(5007) a 


AV .5(5007) 


diameter (in 


Ha) b 




logF(H/3) 


6560? 


WR? 










(A) 




(km/s) 


( A ) 




(km/s) 


50% I max 10% I max 


(cm- 3 ) 


(erg/ s) 






K 5-20 


356.8-03.0 


2007 


Aug 


0.9118 ± 


.0054 


15.09 ± 0.12 


0.639 ± 


.010 


18.02 ± 0.30 


3.2 


6.1 


498 


34.00 


no 


no 


M 1-19 


351.1+04.8 


2005 


May 


0.8781 ± 


.0022 


14.00 ± 0.05 


0.5640 ± 


.0056 


15.56 ± 0.17 


2.6 


5.4 


6370 


34.99 


no 


wels 


M 1-20 


006.1+08.3 


2004 


Jun 


0.7527 ± 


.0037 


9.43 ± 0.08 


0.2969 ± 


.0031 


5.97 + 0.09 


3.0 


4.5 


7500 


35.07 


no 


wels 


M 1-29 


359.1-01.7 


2004 


Jun 


1.0319 ± 


.0033 


18.69 ± 0.08 


0.7693 ± 


.0060 


22.08 ± 0.18 


5.3 


8.8 


3674 


35.34 


yes 


wels 


M 1-31 


006.4+02.0 


2005 


Jul 


0.7954 ± 


.0044 


11.11 ± 0.10 


0.4899 ± 


.0025 


13.11 ± 0.07 


2.2 


4.3 


7830 


35.32 


no 


wels 


M 1-35 


003.9-02.3 


2007 


Aug 


1.0215 ± 


.0063 


18.40 ± 0.14 


0.7336 ± 0. 


.0073 


21.03 ± 0.22 


4.0 


6.7 


6040 


35.10 


no 


no 


M 1-42 


002.7-04.8 


2003 


Jun 


1.0485 ± 


.0038 


19.18 ± 0.09 


0.6772 ± 


.0055 


19.21 ± 0.16 


7.9 


12.4 


1287 


34.95 


yes 


no 


M 1-48 


013.4-03.9 


2005 


Jul 


0.8262 ± 0. 


.0024 


12.23 ± 0.05 


0.5301 ± 0. 


.0046 


14.45 ± 0.14 


3.8 


6.9 


1525 


34.06 


yes 


no 


M 2-13 


011.1 + 11.5 


2006 


Jun 


0.7047 ± 


.0025 


7.25 ± 0.06 


0.4008 ± 


.0020 


10.04 ± 0.06 


2.8 


5.6 


4185 


34.26 


no 


no 


M 2-15 


011.0+06.2 


2006 


Jun 


1.271 ± 


028 


25.23 ± 0.64 


0.937 ± 


.035 


27.3 ± 1.0 


6.6 


9.1 


2130 


34.57 


yes 


no 


M 2-16 


357.4-03.2 


2004 


Jun 


1.1090 + 0. 


.0057 


20.87 ± 0.13 


0.7905 ± 0. 


.0077 


22.74 ± 0.23 


2.5 


4.9 


4625 


35.07 


yes 


no 


M 2-20 


000.4-01.9 


2006 


Jul 


1.2428 ± 


.0041 


24.50 ± 0.09 


0.8262 ± 


.0027 


23.86 ± 0.08 


2.1 


4.3 


2023 


35.02 


no 


WR 


M 2-21 


000.7-02.7 


2005 


May 


1.280 + 0. 


.016 


25.49 ± 0.37 


0.945 ± 


.026 


27.52 ± 0.78 


2.8 


5.3 


10 


34.65 


yes 


wels 


M 2-22 


357.4-04.6 


2007 


Aug 


1.234 ± 


.030 


24.27 ± 0.67 


0.903 ± 


.039 


26.3 ± 1.1 


3.1 


6.9 


2300 


34.18 


yes 


no 


M 2-23 


002.2-02.7 


2004 


Jun 


0.855 ± 0.0021 


13.22 ± 0.05 


0.5258 ± 


.0030 


14.30 ± 0.09 


2.5 


4.3 


10600 


35.00 


no 


no 


M 2-26 


003.6-02.3 


2006 


Jul 


1.319 ± 0. 


.024 


26.50 ± 0.56 


1.011 ± 


.033 


29.55 ± 0.99 


8.2 


11.2 


452 


34.14 


yes 


no 


M 2-27 


359.9-04.5 


2004 


Jun 


0.9625 ± 


.0034 


16.64 ± 0.08 


0.6230 ± 


.0044 


17.46 ± 0.13 


1.9 


6.4 


7360 


35.11 


no 


wels 


M 2-29 


004.0-03.0 


2004 


Jun 


0.7864 + 


.0026 


10.77 + 0.06 


0.4771 ± 


.0036 


12.68 ± 0.11 


3.6 


7.1 


2670 


34.53 


no 


no 


M 2-30 


003.7-04.6 


2004 


Jun 


1.184 ± 


.011 


22.91 ± 0.25 


0.864 ± 


.016 


25.02 + 0.47 


3.2 


5.6 


2268 


34.83 


yes 


wels 


M 2-31 


006.0-03.6 


2004 


Jun 


1.1907 ± 


0068 


23.10 ± 0.16 


0.8266 ± 


0060 


23.87 ± 0.18 


2.3 


5.6 


4710 


35.19 


no 


WR 


M 2-33 


002.0-06.2 


2007 


Aug 


0.6988 ± 0. 


.0031 


6.98 ± 0.07 


0.3441 ± 


.0046 


7.98 ± 0.14 


4.1 


6.1 


1259 


34.94 


no 


wels 


M 2-39 


008.1-04.7 


2006 


Jul 


1.0341 ± 


.0053 


18.76 ± 0.12 


0.5258 ± 


.0063 


14.31 ± 0.19 


2.4 


3.6 


3165 


34.62 


no 


wels 


M 2-4 


349.8+04.4 


2007 


Aug 


0.8088 ± 


.0041 


11.63 ± 0.09 


0.5025 ± 0. 


.0035 


13.58 ± 0.10 


1.7 


4.4 


6970 


35.39 


no 


no 


M 2-8 


352.1+05.1 


2006 


Jun 


0.9742 ± 


.0019 


17.00 ± 0.04 


0.6275 ± 


.0021 


17.61 ± 0.06 


2.2 


5.3 


1670 


34.44 


yes 


WR 


M 3-10 


358.2+03.6 


2004 


Jun 


1.0963 ± 


.0085 


20.52 ± 0.19 


0.7351 ± 


.0069 


21.01 ± 0.21 


2.5 


4.7 


11497 


34.99 


yes 


no 


M 3-14 


355.4-02.4 


2004 


Jun 


1.1585 ± 


.0075 


22.23 ± 0.17 


0.8005 ± 


.0087 


23.06 ± 0.26 


3.5 


7.5 


3180 


35.28 


yes 


no 


M 3-15 


006.8+04.1 


2004 


Jun 


0.9783 ± 


0033 


17.12 ± 0.08 


0.6470 ± 


0062 


18.23 ± 0.19 


4.0 


6.1 


9450 


35.67 


no 


WR 


M 3-16 


359.1-02.3 


2005 


May 


1.0279 ± 


.0074 


18.58 ± 0.17 


0.6229 ± 


.0085 


17.46 ± 0.25 


6.0 


11.2 


1155 


35.06 


no 


no 


M 3-20 


002.1-02.2 


2007 


Aug 


1.1206 ± 0. 


.0093 


21.20 ± 0.21 


0.7629 ± 


.0046 


21.89 ± 0.14 


3.5 


4.8 


10 


34.64 


no 


wels 


M 3-21 


355.1-06.9 


2004 


Jun 


0.8936 ± 


.0062 


14.49 ± 0.14 


0.5682 ± 


.0056 


15.69 ± 0.17 


1.8 


4.4 


3172 


35.06 


no 


no 


M 3-26 e 


004.8-05.0 


2005 


Sop 


1.462 ± 


.052 


30.2 ± 1.2 


1.141 




33.5 


7.0 


10.7 


1415 


34.13 


yes 


no 


M 3-32 


009.4-09.8 


2005 


Jul 


1.176 ± 


.012 


22.70 ± 0.27 


0.811 ± 


.013 


23.37 ± 0.38 


4.0 


7.4 


1850 


34.38 


yes 


no 


M 3-33 


009.6-10.6 


2004 


Jun 


1.098 ± 


.012 


20.55 ± 0.28 


0.806 ± 


.020 


23.24 ± 0.61 


5.3 


8.8 


1140 


34.05 


yes 


wels 


M 3-38 


356.9+04.4 


2004 


Jun 


0.9558 ± 


.0086 


16.44 ± 0.20 


0.5289 ± 


.0077 


14.41 ± 0.23 


3.2 


4.4 


9550 


34.31 


yes 


no 


M 3-42 


357.5+03.2 


2003 


Jun 


1.710 ± 


.096 


36.4 ± 2.2 


0.602 ± 


.046 


16.7 ± 1.4 


4.1 


9.8 


1000 


33.72 


yes 


no 


M 3-45 


359.7-01.8 


2005 


Jul 


1.1658 ± 0. 


.0068 


22.43 ± 0.16 


0.8117 ± 


.0095 


23.41 ± 0.28 


5.1 


7.4 


7897 


34.86 


yes 


no 


M 3-54 


018.6-02.2 


2006 


Jul 


1.340 ± 


.020 


27.03 ± 0.46 


0.957 ± 


.026 


27.89 ± 0.79 


4.4 


7.4 


1385 


33.97 


yes 


no 


M 4-3 


357.2+07.4 


2005 


May 


0.8502 ± 0. 


.0022 


13.07 + 0.05 


0.5271 ± 


.0025 


14.35 ± 0.07 


3.6 


6.0 




35.27 


no 


no 


M 4-6 


358.6+01.8 


2004 


Jun 


1.1311 ± 


.0067 


21.48 ± 0.15 


0.772 ± 0. 


.014 


22.15 ± 0.40 


3.7 


6.6 


12881 


34.26 


no 


no 


M 4-7 


358.5-02.5 


2006 


Jun 


1.0558 ± 0. 


.0060 


19.38 ± 0.14 


0.730 + 0. 


.012 


20.86 ± 0.36 


4.8 


8.0 


1580 


34.15 


no 


no 


PC 12 


000.1 + 17.2 


2005 


May 


0.8346 ± 


.0019 


12.53 ± 0.04 


0.4749 ± 


.0019 


12.61 ± 0.06 


3.1 


6.3 


5110 


34.27 


no 


no 


Te 1580 


002.6+02.1 


2007 


Aug 


1.883 ± 


.065 


40.6 ± 1.5 


1.367 ± 


.075 


40.42 ± 2.23 


7.4 


9.9 


1200 


34.48 


yes 


no 



a This is the observed line width, uncorrected for instrumental, thermal, or fine structure broadening. 



These diameters were measured by collapsing the Ha spectra into one-dimensional spatial profiles, and measuring the diameters at 50% and 10% of peak intensity. 
c Thc electron densities arc measured from the [S II] AA6716, 6731 lines. When no density is given, this information is unavailable. 

d Classifications of "WR" and "wels" are taken from an extended version of Table 4 from[G6rny ct al. (20Q4j) kindly provided by S. Gorny. A classification of "no" indicates that no 
information is available or that the central star is classified as non-WR. 

e The [O III]A5007 line width was measured using IRAF's splot. 



CO 



Table 2 
Statistical Tests (U-test) 



Hypothesis 


Prob(Fals 


*0 


Comment 


AV .5,5007(hottcr CS) > AV .5,5007(cooler CS) 


1.3 X 10" 


7 


H 1-67 excluded 


AV . 5, He (hotter CS) > AVa. 5, Ha (cooler CS) 


7.4 x 10" 


10 




diamcter(hotter CS) > diameter(cooler CS) 


7.0 x 10" 


4 


measured at 10% intensity 


density (cooler CS) > density (hotter CS) 


8.0 x 10" 


3 


based upon[S n]AA6716, 6731 


L(H/3, cooler CS) > L(H/3, hotter CS) 


2.2 x 10" 


■4 
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Spectral Line 1 = 6560. 1i 




6560 6562 6564- 6566 6568 6570 



-wHe II 6560 absent 
-■-He II 6560 present 














/K A 













10000 15000 20000 

density (cin-3) 



Fig. 2. — The Ha line profile for M 2-16 is an 
example of our spectra. The crosses are the data 
and the solid line is the fit. The faint line on the 
blue wing of Ha is the He II A6560 line and has an 
intensity of 0.9% that of Ha. For the Ha line, the 
S/N in the flux is 200, which is approximately the 
median value for our sample. 



Fig. 5. — The distributions of densities are very 
different for the Bulge planetary nebulae hosting 
cooler and hotter central stars (He IIA6560 ab- 
sent and present, respectively). The probability 
of drawing these distributions of densities at ran- 
dom from the same parent population is 8.0 x 10~ 3 
(Table (2). 
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diameter (10% I m .„) 

Fig. 4. — The distributions of diameters at 10% of 
peak intensity (or 50% of peak) are very different 
for the Bulge planetary nebulae divided accord- 
ing to whether the He IIA6560 line is absent or 
present. The probability of drawing such differ- 
ent distributions of diameters at random from the 
same parent population is 7.0 x 10 -4 (Table [2]). 




logL(Hf)) (erg/s) 

Fig. 6. — As was found for the line widths, diam- 
eters, and densities, the distributions of H(3 lumi- 
nosities are very different for the Bulge planetary 
nebulae hosting cooler and hotter central stars 
(He IIA6560 absent and present, respectively). 
The probability of drawing these distributions of 
luminosities at random from the same parent pop- 
ulation is 2.2 x 10~ 4 (Table©. 
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Fig. 1 . — The selection criteria for our sample of Bulge planetary nebulae were designed to yield a sample that 
mimics the populations of bright extragalactic planetary nebulae in bulge-like systems. The extragalactic 
planetary nebulae arc shown in solid symbols while the Bulge planetary nebulae, divided according to the 
presence or absence of He IIA6560 are shown with open symbols. The dotted lines indicate the limits we chose 
for the H/3 luminosity and [O IH]A5007/H/3 ratio were motivated by the range of these parameters among 
extragalactic planetary nebulae with spectroscopic observations dJacobv fc Ciardullo 1ll999t iRicher et al 



1999; Walsh ct al. 1999; Dudziak ct al. 2000; R oth et al. 1 12004: Mcn dez et al. 2005; Zi ilstraa et al. 200' 



Goncalves et al. 1 120071: IRicher fc McCall 1120081) 



0G; 
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Fig. 3. — We find distributions of line widths that are very different if we divide our objects according to 
whether the He IIA6560 line is absent or present. The distribution of Ha line widths is on the left while the 
distribution of [O III]A5007 line widths is on the right. The probability of drawing such different distributions 
(He IIA6560 absent or present) at random from the same parent population is less than 1.3 x 10" 7 (Table 
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